Detailed abundances of Red Giants in the Globular Cluster 
NGC 1851: C+N+O and the Origin of Multiple Populations 



S. Villanova and D. Geisler 

Departamento de Astronomia, Casilla 160, Universidad de Concepcion, Chile 
svillanova@astro-udec . cl 
and 

G. Piotto 

Dipartimento di Astronomia, Universitd di Padova, Vicolo dell'Osservatorio 3, 1-35122 Padua, Italy 



ABSTRACT 

We present chemical abundance analysis of a sample of 15 red giant branch (RGB) stars 
of the Globular Cluster NGC 1851 distributed along the two RGBs of the (v, v-y) CMD. We 
determined abundances for C+N+O, Na, a, iron-peak, and s-elements. We found that the two 
RGB populations significantly differ in their light (N,0,Na) and s-element content. On the other 
hand, they do not show any significant difference in their a and iron-peak element content. 
More importantly, the two RGB populations do not show any significant difference in their total 
C+N+O content. 

Our results do not support previous hypotheses suggesting that the origin of the two RGBs 
and the two subgiant branches of the cluster is related to a different content of either a (including 
Ca) or iron-peak elements, or C+N+O abundance, due to a second generation polluted by SNell. 

Subject headings: globular clusters: general — globular clusters: individual(NGC 1851) 



1. Introduction 

Recent investigations have fostered new inter- 
est in globular cluster (GC) stellar populations. 
The traditional vision of GCs hosting simple, sin- 
gle age, single metallicity stars has been shattered 
by the discovery of multiple sequences in the color- 
magnitude diagrams (CMD) and of large disper- 
sions in light or heavy element abundances, often 
showing well defined patterns, which are the re- 
sult of nuclear processes in a previous generation 
of stars. 

NGC 1851 is one of the GCs where multiple stel- 
lar populations have been clearly identified. In- 
deed, in recent years, this GC has been the sub- 
ject of a wide photometric and spe ctroscopic ob 
servat ional campaign. First of all, iMilone et al. 



presence o f at least two stellar g enerations. As dis- 
cussed by IMilone et al. I (|2008[) . if the SGB split 
would be due to a different age between the two 
populations with similar metal content, the sec- 
ond generation of stars should have formed ~ 1 
Gyr after the first one. However other factor s 
could be at work. As pointed out bv lCassisi et ah 



(2008) have identified a double subgiant branch 
(SGB) in the cluster CMD, clearly suggesting the 



1)20081 ) and lVentura et al. I (|2009h . the SGB split- 
ting could also be explained by a difference in 
the total C+N+O content, with an abundance for 
the fainter SGB 2-3 times larger than that of the 
brighter one. In this case, the age spread would 
be smaller, of the order of a few hundred Myrs. 

Unfortunately, we have scarce information 
on detailed abun dances of NGC 1851 stars. 
Hesser at al. I (Il982h collected low resolution spec- 
tra for 18 stars. They found a constant iron con- 
tent, but 3 out of 8 of their bright RGB stars 
show anomalously strong CN bands. More re- 
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cently, lYong fc Grundahl I ||2008j) analyzed UVES 
spectra of eight giants. With all the limita- 
tions coming from the s mall number statistics, 
Yong fc Grundahfl ( 20081) results seem to indicate 
the presence of a Na-0 anticorrelation, and the 
presence of a large star to star abundance scatter 
of s-process elements Zr and La, with some hints 
of a bimodal distribution. They also find a large 
star to star variation in the strength of the Ba 
lines. Th ey also show tha t the s plit of the RGB 
found by ICalamida et~al~l (j2007l ) in a Stromgren 



photometry CMD is related to the double SGB. 
Interestingly, the fraction of La/Zr-strong stars 
(40%) is tantalizing similar to the fraction of SGB- 
faint stars (45%), and to the fraction of the CN- 
st rong, Ba-strong, and Sr-strong RGB stars (40%) 



of iHesser at al.l (|1982T) . But the sample statistics 
for stars with abundance esti mates is very poor. 
Recently, lYong et al. I (l2009h suggested the pres- 
ence of a large C+N+O spread (~0.6 dex) among 
NGC 1851 giant stars. Their work is based on 
high resolution spectra of four stars. This spread 
could imply that one population has 4 times larger 
C+N+O content with respect to t he other, sup- 



portin g th e scenario prop o sed b y ICassisi et al 



(|2008l) and I Ventura et al. I (|2009h . However, also 



in this case, the statistics is very poor. 

In order to account for such observational re- 
sults, it has been proposed that a second gen- 
eration of stars may have formed from a resid- 
ual reservoir of gas polluted by intermediate 
mass asymptotic giant branch (AGB) star ejecta 
(jP'Antona et al. Il2002l ) , or from the ejecta of m as- 
sive fast rotating stars dDecressin et al. I 2007 ). 
A very recent work by iLee et al. I (|2009j ) proposes 
another possibility. Using Stromgren photome- 
try and the Ca narrow-band filter cente r ed on 
the H-K Ca lines at -3950 A, ILee et al. I (j2009l ) 



find a significant color split in the RGB of many 
clusters, incl uding NGC 1851. Because of the ob- 



served split, ILee et al. I (2009) propose that the 
two RGBs have a different Ca content, of the 
order of 0.3 dex. In this case, the second genera- 
tion of stars could originate from gas polluted by 
SNell explosions, the only mechanism able to pro- 
duce a Ca enhancement. However, this scenario 
is quite controversial, because multiple CMD se- 
quences have been found in other clusters, where 
abundances based on high resolution spectra give 
a negligible or very low Ca spread. For exam- 
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Fig. 1. — V vs. V-I and v vs. v-y CMDs of our 
target stars. The double RGB is visible in the v 
vs. v-y CMD but not in the V vs. V-I CMD. The 
blue RGB is indicated as a blue line, the red RGB 
as a red line. Ba-poor stars are plotted as filled 
circles, while Ba-rich as open circles. 



pie, th ere i s the case of M4 , where IMarino et al 



(|2008l ) and ICarretta et al. 
spread at all, while ILee et al. 



( 2010aD fo und no Ca 
20091 ) presented 



a CMD with a double RGB. In M4, a clear RGB 
spread is visible in ultraviolet broad-band photom- 
etry (see Fi g 11 of IMarino et al~l (120081 )). but, as 
discussed in Marino et al. ( 20081 ). the photomet- 
ric spread can be better explaind by a spread in 
C,N,0. 

The only way to solve the enigma is to mea- 
sure C,N,0, a, iron-peak, and s-process element 
abundances for a statistically significant sample of 
RGB stars. For this purpose we observed a sam- 
ple of 15 RGB stars distributed among the two 
RGBs of the Stromgren (v,v-y) CMD with high 
and medium resolution spectroscopy. The obser- 
vations, data reduction, chemical abundance anal- 
ysis, and results will be presented in the following 
sections. 

In this letter we discuss only the results on Fe, 
Ca, s-process, and C+N+O content, i.e. the most 
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relevant elements in the context of the present de- 
bate on multiple stellar populations in GCs. Full 
details will be presented in a subsequent paper. 
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Fig. 2.— [Ba/Fe] vs. [Fe/H] and [Y/Fe] vs. 
[Zr/Fe] for our target stars. Ba-poor stars are 
plotted as filled circles, while Ba-rich as open cir- 
cles. Red crosses give mean values and errors of 
the mean for the two groups of stars. 



2. Observations 

Our dataset consists of high and medium res- 
olution spectra collected in November-December 
2008. The spectra come from single ~1500 and 
—500 s exposures, obtained with the MIKE and 
Mage slit spectrographs respectively, mounted at 
the Magellan/Clay telescope. Weather conditions 
were excellent with a typical seeing of 0.5 arcsec. 
We selected 15 isolated stars from 0.5 magnitudes 
above the RGB-bump to the tip of the RGB, in 
the magnitude range 15.5<F<13.5. The MIKE 
spectra have a spectral coverage from 4900 to 8700 
A with a resolution of - 32000 and S/N>100 at 
6300 A. The Mage spectra have a spectral cover- 
age from 3700 to 9500 A with a resolution of ~ 
5000 and S/N>50 at 4300 A. 



Data were reduced using IRAF@ including bias 
subtraction, flat-field correction, wavelength cali- 
bration, scattered-light and sky subtraction, and 
spectral rectification. 

In the case of MIKE observations, spectra were 
corrected for fringing using the spectrum of a hot 
fast rotating star, which gave the best result com- 
pared with other methods. However, fringing did 
not affect the spectral region we were interested in 
(below 8100 A). 

Targets were carefully selected in order to have 
stars in each of the two distinct RGBs clearly sep- 
arated in the high quality v vs. v-y Stromgren 
photometry CMD (see Fig. [TJ lower panel) kindly 
provided by F. Grundahl. In the following, we 
used also VI Cousins photometry from Y. Mo- 
many. Note that there is no RGB split in the V vs. 
V-I CMD (Fig. HJ upper panel). Filled and open 
circles in Fig. [T]are Ba-poor and Ba-rich stars, re- 
spectively (as measured from our spectra, see Sec. 
4). All stars were observed with both spectro- 
graphs. 

Radial velocities were measured by the fxcor 
package in IRAF, using a synthetic spectrum as 
a template. The mean heliocentric value for our 
targets is 320.3± 1.1 km/ s , whi le the dispersion 
is 4.4±0.8 km/s. lHarris I (jl996l) gives 320.5±0.6 
km/s as heliocentric radial velocity for NGC1851, 
and the typical dispersion for a cluster of its mass 
is -4-5 km/s (jPrvor fc Mevlan I Il993h . All our 
targets have a radial velocity within ±15 km/s of 
the cluster radial velocity. On the basis of this re- 
sult, we conclude that all our targets are cluster 
members. 

3. Abundance analysis 

The chemical abundances for Na, Si, Ca, Ti, Fe, 
and Ni were obtained from the eq uivalent widths 
EWs) of the spectral lines 



See Marino et al. 



20081 ) for a more detailed explanation of the 
method we used to measure the EWs and for the 
adopted solar abundances. For the other elements 
(C, N, O, Y, Zr, Ba), whose lines are affected by 
blending, we used the spectral synthesis method. 
Na presents few features in the spectrum, so in 



1 Sec www.lco.cl 



2 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation 
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this case abundances derived from the EWs were 
cross-checked with the spectral synthesis method 
in order to obtain more accurate measurements. 
MIKE spectra were used to obtain abundances for 
all the elements but C, which was obtained from 
Mage data. Only lines not contaminated by tel- 
luric lines were used. 
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Fig. 3. — Ni and a-element abundances of our 
target stars. Ba-poor stars are plotted as filled 
circles, while Ba-rich as open circles. Red crosses 
are mean values and errors of the mean for the two 
groups of stars. 

Atmospheric parameters were obtained in the 
following way. First of all T e ff wa s derived from 



the V -I color using the relation b y lAlonso et al 



( 19991 ) and the reddening from Harris (1996) 



Surface gravities (log(g)) were obtained from the 
canonical equation: 



log(^) = log(^)+4.1og& 



log(^) 



where the mass M/Mq was derived from 



isochrone fitting using the Padova database (jMarigo et al. 
2008), and the luminosity L/Lq was obtained 
from the absolute magnitude My assuming an 
apparent dist ance modulus of (m — M )v=15.47 
(jHarris I fl996). The bolometric correction (BC) 
was derived by adopting the relation BC-T c ff from 



1999). Finally, microturbulence 
obt ained from the relation of 
Gratton et al. 1 (|l996l ). 



Alonso et al 
velocity (vt, 



was 



These atmospheric parameters were considered 
as initial estimates and were refined during the 
abundance analysis. As a first step at mospheric 
models were calculated using ATLAS9 (jKurucz I 
1970f ) and assuming the initial estimate of T e g, 
logfq) , and vt, and the [Fe/H] value from lHarris I 



fll996| . 

Then T c fj, w t , and log(g) were adjusted and new 
atmospheric models calculated in an interactive 
way in order to remove trends in Excitation Po- 
tential (E.P.) and equivalent widths vs. abun- 
dance for T e g and v t respectively, and to satisfy 
the ionization equilibrium for log(g). Fel and Fell 
were used for this purpose. The [Fe/H] value of 
the model was changed at each iteration according 
to the output of the abundance analysis. 
The Local Ther modynamic Eq uilibrium (LTE) 
program MOOG (jSneden Ill973h was used for the 
abundance analysis. 

We checked the reliability of our atmospheric 
parameters in the following way. Intrinsic V-I 
colours of our st ars were obtained fro m our T e g- 
by inverting the Alonso et al. ( 19991 ) equation. 



Comparing intrinsic and measured colours we got 
a reddening of E(B-V)=0.01±0.01 (assuming E(V- 
I)=1.24xE(B-V)). 

This y alue is in very good agreement with lHarris I 
(|l996l ). that gives E(B-V)=0.02. 
A further check to test the reliability of our T c g 
and log(g) values was performed on Arcturus, 
which is an important reference for every study 
on RGB stars. 

In order to have statistically independent estima- 
tions, we measured the atmospheric parameters on 
3 high- resolution high-S/N spectra obtained with 
HARPS, the high precision spectrograph mounted 
at the 3.6m telescope in La Silla. 
We obtained T off =4290±ll K and log(g)= 
2.00±0.02 as mean values and internal errors, 

which agree extremely well with li terature values 

(T Rff =4290±30 K, log(gWl.9±0.1. lGriffin fc Lvnas-Grav 



1999). 



Our conclusion is that our T c g- and log(g) scales 
appear free from significant systematic errors. 

The linelists for the che mical analysis were ob - 
tained from many sources (jGratton et al. I (|2003l ) , 
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McWilliam 



VALD & NIST^ jMcWilliam et al. I (Il994h 
(|l998l ). SPECTRUJV0, and SCAnI), and cali 



brated using the Solar-inverse technique by the 
spectral synthesis method. For this purpose we 
used the hig h resolution, high S/ N NOAO So- 



lar spectrum ([Kurucz et al. 1 119841 ) . For Ba lines 
we took the hyperfine splitting into account. We 
emphasize the fact that all the linelists were cal- 
ibrated on the Sun, including those used for the 
spectral synthesis. 

In particular, our determinations of C,N,0 abun- 
dances are based on the G-band at 4310 A, the 
CN lines at 8003 A, and the forbidden O line at 3 
6300 A respectively. CN lines at 8003 A were used 8 
also to estimate the C 12 /C 13 ratio. J 

For the Sun we have loge(C)=8.49, loge (N)=7.95, ~ 

and lo ge(Q)=8.83, values that agree with lGrevesse fc Sauval 
(1998) within 0.03 dex. 





These three features were also checked on the 
hig h resolution, h igh S /N spectrum of Arcturus 
bv iHinkle et al. (2003) in order to verify if fur- 
ther blends affecting cool RGB stars but not the 
Sun were present. We found no significant con- 
tamination. 

We obtained the following results for Arcturus: 
[C/Fe]=-0.10±0.05, [N/H]=+0.34± 
0.05, [Q/Fe1=+0.38±0.0 5 (internal errors ). Thes e 
values well agree with IPeterson et al. I (1993), 
which gives [C/Fe]=+0.0±0.1, [N/H]=+0.3±0.1, 
[O/Fc]=+0.4±0.1. We consider this as a test 
which proves the reliability of our linelists for 
C,N,0 not only on Sun-like stars, but also on cool 
RGB stars. 

Abundances for C, N, and O were determined all 
together in an interactive way, because for the 
temperature of our stars, carbon, nitrogen, and 
oxygen form molecules and as a consequence their 
abundances are related to each other. 
Our C 12 /C 13 estimation (C 12 /C 13 <15 for all our 
targets) confirm that the stars are affected by evo- 
lutionary mixing, as expected from their position 
in the CMD. This can affect the primordial C,N,0 
abundances separately, but not the total C+N+O 
content because these elements are transformed 
one into the other during the CNO cycle. 



0.2 0.4 
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3 See |http://vald.astro.univie.ac.at/~vai d/php/vald.php 
and http://physics.nist.gov/PhysRcfData/ASD/lincs_form.html 

4 See |http:// www.phys.appstate.edu/spectrum/spectrum.html 
and references therein 

5 See |http://www. astro.ku.dk/~ uffegj / 



Fig. 4. — Na-0 anticorrelation and CNO abun- 
dances for our targets. Ba-poor stars are plotted 
as filled circles, while Ba-rich as open circles. Red 
crosses are mean values and errors of the mean for 
the two groups of stars. 



4. Results and Discussion 

First of all, in Fig. [2] we plot the abundance of 
Ba vs. Fe. (upper panel). In this plot we clearly 
see that the Ba distribution is bimodal, while the 
iron abundance for the two Ba groups is the same 
within the errors. In the following analysis, we 
divided our stars into two groups, one Ba-poor 
(filled circles in the figures) , and the other Ba-rich 
(open circles). 

The individual and mean abundances we obtained 
for the two groups are summarized in Tab. Q] For 
the mean abundances we reported also the internal 
errors (standard deviation of the mean). 

The two groups have: 

[Fc/H]=-1.23±0.01, [Ba/Fe] = +0.08±0.02 
and 

[Fe/H]=-1.22±0.01, [Ba/Fc]=+0.51±0.03 

as mean values respectively. The lower panel 
of Fig. [2] displays the abundances of two other s- 
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Table 1 

Abundances of the two groups of stars. Mean abundaces of the two populations and the 
related internal errors are reported at the bottom of the table. errors associated 

WITH THE MEASUREMENTS ARE PLOTTED IN FlGS. I2ll3l AND |H 



Id 


[Fe/H] 


[Ba/Fe] 


[Zr/Fe] 


[Y/Fe] 


[Ni/Fe] 


[Si/Fe] 


[Ca/Fe] 


[Ti/Fe] 


[Na/Fe] 


[O/Fe] 


C±N±0 


Ba-poor stars 


13 


-1.24 


0.06 


0.16 


-0.01 


-0.10 


0.23 


0.33 


0.27 


-0.26 


0.26 


7.92 


30 


-1.23 


0.06 


0.05 


0.17 


-0.12 


0.31 


0.32 


0.26 


0.06 


0.18 


8.09 


31 


-1.26 


-0.02 


0.11 


0.11 


-0.08 


0.23 


0.37 


0.30 


-0.31 


0.37 


8.00 


43 


-1.22 


0.06 


0.14 


0.17 


-0.11 


0.39 


0.31 




0.34 


-0.13 


7.98 


51 


-1.19 


0.14 


0.22 


0.12 


-0.07 


0.30 


0.32 


0.30 


0.26 


-0.03 


8.06 


53 


-1.22 


0.13 


0.15 


0.09 


-0.10 


0.32 


0.42 


0.28 


-0.28 


0.23 


7.92 


G3 


-1.26 


0.14 


0.12 


0.14 


-0.03 


0.35 


0.39 


0.31 


0.18 


-0.07 


7.90 


68 


-1.24 


0.13 


0.13 


0.16 


-0.07 


0.25 


0.35 


0.26 


0.36 


-0.11 


8.01 


Ba-rich stars 


35 


-1.24 


0.48 


0.26 


0.23 


-0.08 


0.38 


0.37 


0.28 


0.31 


-0.09 


7.80 


14 


-1.21 


0.54 


0.32 


0.22 


-0.05 


0.32 


0.38 


0.26 


0.20 


0.10 


8.01 


1G 


-1.24 


0.50 


0.16 


0.21 


-0.08 


0.41 


0.41 


0.28 


0.40 


-0.09 


8.03 


18 


-1.19 


0.47 


0.25 


0.25 


-0.06 


0.28 


0.39 


0.25 


0.73 


-0.49 


8.12 


20 


-1.24 


0.41 


0.23 


0.23 


-0.07 




0.42 


0.31 


0.57 


-0.19 


8.08 


8 


-1.18 


0.59 


0.27 


0.22 


-0.11 


0.34 


0.32 


0.27 


0.44 


-0.22 


8.11 


9 


-1.22 


0.63 






-0.05 


0.38 


0.37 




0.61 


-0.33 


7.97 


Mean abundances and internal errors 


Ba-poor 


-1.23 


+0.09 


±0.14 


±0.12 


-0.09 


±0.30 


±0.35 


±0.28 


±0.04 


±0.09 


7.99 




±0.01 


±0.02 


±0.02 


±0.02 


±0.01 


±0.02 


±0.01 


±0.01 


±0.10 


±0.07 


±0.02 


Ba-rich 


-1.22 


±0.52 


±0.25 


±0.23 


-0.07 


±0.35 


±0.38 


±0.28 


±0.47 


-0.19 


8.02 




±0.01 


±0.03 


±0.02 


±0.01 


±0.01 


±0.02 


±0.01 


±0.01 


±0.07 


±0.07 


±0.04 
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elements: Y and Zr. Again the two groups of Ba- 
rich and Ba-poor stars are separated, in the sense 
that the Y and Zr content is significantly different 
for Ba-poor and Ba-rich stars. The two groups 
have [Zr/Fe]=+0.14±0.02, [Y/Fe]=+0. 12+0.02 
and [Zr/Fe]=+0.25±0.02, [Y/Fe]=+0.23±0.01 re- 
spectively. The results displayed in the lower panel 
of Fig.[2]can be considered as an independent con- 
firmation of the presence of two groups of stars 
with different s-element content in NGC 1851. 

An interesting result comes from Fig. [1] (lower 
panel) which clearly shows that all Ba-poor stars 
are located in the bluer RGB, while all Ba-rich 
stars are located in the redder one. The two se- 
quences seem to be well separ ated. As in the case 



age [Na/Fe]=+0.04±0.10and [O/Fc]=+0.09±0.07, 
while Ba-rich stars have [Na/Fe]=+0.47±0.07, 
[O/Fe]=-0.19±0.07. 

We found that Na is correlated also with C and 
N. The two groups of stars have different mean 
N contents ([N/Fe] = 1.09±0.08 dex for Ba-poor 
stars, [N/Fe]=0.73±0.16 for Ba-rich stars, but 
very similar mean C content ([C/Fe]=-0.64±0.09 
dex for Ba-poor stars, [C/Fe]=-0.72±0.08 for Ba- 
rich stars). However, as in the case of M22, C, N 
and O do not show a bimodal distribution (Vil- 
lanova et al., in preparation). How C,N,0 content 
is related to the double RGB visible in the (v, v-y) 
CMD is not clear. We know that the v band con- 
tai ns many CH and CN features which, according 



of M22 (|Marino et al. II2009I) . our result implies to |Marino et al. | (|2008|), can contribute to split 

that s-elements allow us define two distinct popu- thc RGB in UV or blue CMDs. CH b ands affect 

lations also in NGC1851. also the Ca narrow-band filter used by Lee et al. 

and this could explain their RGB split 
without invoking a Ca spread. 
Any interpretation of the CMD and of the spec- 
troscopic results involving C,N,0 should also take 
into account the evolutionary stage of our targets. 
The investigation of light-element correlations and 
the relation between the double RGB and the 
C,N,0 content is beyond the scope of this letter 
and will be discussed in detail in a future paper. 

Our final interesting result concerns C+N+O. 
Fig. |4] (lower panel) shows that the two groups 
of stars have the same C+N+O content, within 
la. In particular, the Ba-poor group have 
loge(CNO)=7.99±0.02, while the Ba-rich have 
loge(CNO)=8. 02+0.04. 

About the mean C12/C13 content of the two pop- 
ulations we found they have the same values (^5) 
within the errors. We expected the Ba-rich stars 
to have a lower C12/C13 value as the results of 
the contamination of the pristine gas they formed 
from by highly processed ejecta. 
However we point out that our results about 
C12/C13 are preliminar. Also the evolutionary 
stage of the targets must be taken into account, 
which could alter the primordial mean C12/C13 
ratio of the two populations making them to as- 
sume a similar value. 

A better estimation of C12/C13 for each targets 
and a theoretical treatment of the expected differ- 
ence between the two population will be presented 
in a future paper. 

Figure[4]excludes the possibility that the second 



Now, it is instructive to analyze the a-element 
content of the two groups. We included in this 
analysis also Ni, another iron-peak element. For 
this purpose, we plot in Fig. [3] Ni vs. Ba (upper- 
left panel), Si vs. Ba (upper-right panel), Ca vs. 
Ba (lower- left panel), and Ti vs. Ba (lower-right 
panel) . 

Thc average value of Ni, Si, Ca, and Ti for the two 
different Ba groups is summarized in Table [TJ with 
the corresponding internal dispersion of the mean 
plotted in the figure. As already found for iron 
(Fig. [5]), the Ba-rich and Ba-poor groups have no 
appreciable difference in their iron-peak element 
and a-element content, which have similar values 
within l-1.5cr. 

This result is in contra st with the hypothe- 
. (l2009h . and consistent with 
"(l2010alT s statement that most 



sis by iLee et al. 



Carretta et al. 



GCs do not have an appreciable spread in their 
iron-peak or a-elcmcnt content. Our results also 
imply that in the case of NGC1851 pollution by 
SNell cannot explain the double population. 

As for the lighter elements, Fig. [4] (upper panel) 
confirms that the cluster has a very well de- 
fined Na-0 anticorrelatio n, as already shown by 
Yong fc Grundahll (l2008h . 

However, here we have an interesting new re- 
sult. The two Ba groups are partially separated 
in the Na vs. O plot, with some overlap at 
—0.1 <[0/Fe]< +0.1. Ba-poor stars have an aver- 
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generation of stars have a significa ntly enhanced 
C+N+O content, as suggeste d by ICassisi et al 



2008|), IVentura et al~l (2009t ). and lYong et al. 



20091 ) . U nfortunatly w e do no t have stars in com- 



mon with Yong et al. (l2009f) for a comparison 



Here we note that lYong et al. ] (120091) targets are 
very bright, close to the RGB-tip. It is well possi- 
ble that some of their stars are indeed AGB stars. 
In such a case, the chemical pattern shown by 
those stars may have been highly altered by in- 
ternal nuclear processes. 
However, as suggested by lCassisi et al. 



(2008|), if 



the first generation had loge(CNO)=8.00 dex, then 
a second generation with double C+N+O content 
would correspond to loge(CNO)=8.30 dex. At odd 
with that, a difference in loge(CNO) of 0.03, as 
we find, implies that the Ba-rich population have 
about 1.1 times the C+N+O content of the Ba- 
poor one, which is consistent with a null difference 
when we consider the measurement errors. 
The dicothomy in Ba and other s-process element 
content, as well as the clear spread in the CMD 
of Fig. [1] clearly points towards the presence of 
two distinct st ellar populations in N GC 1851, as 



found also by iMilone et al~l (J2OO8) in the SGB 



However, the constancy of the CNO content casts 
doubts on the scenario of a second population born 
from material polluted by first generation inter- 
mediate mass stars. The possibility of a merge of 
two distinct clusters, as very recently suggested by 
Carretta et al~l (|2010bl) should be furher explored. 
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